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lattanziOto . astro . it 
ABSTRACT 

Modern astrometry is based on angular measurements at the micro- arcsecond 
level. At this accuracy a fully general relativistic treatment of the data reduction 
is required. This paper concludes a series of articles dedicated to the problem 
of relativistic light propagation, presenting the final microarcsecond version of a 
relativistic astrometric model which enable us to trace back the light path to its 
emitting source throughout the non-stationary gravity field of the moving bodies 
in the Solar System. The previous model is used as test-bed for numerical com- 
parisons to the present one. Here we also test different versions of the computer 
code implementing the model at different levels of complexity to start exploring 
the best trade-off between numerical efficiency and the //as accuracy needed to 
be reached. 

Subject headings: astrometry — gravitation — reference systems — relativity — 
time 



1. Introduction 

Modern space technology will soon provide stellar positioning with micro-arcsecond ac- 
curacy (//as). At this level one has to take into account the general relativistic effects on light 
propagation arising from metric perturbations due not only to the bulk mass but also to the 
rotational and translational motion of the bodies of the Solar System, and to their multipole 
structure (see Kopeikin and Mashhoon 2002; Klioner 2003; Le Poncin-Lafitte et al. 2004 
and references therein). Our aim is to develop a Relativistic Astrometric MODel (RAMOD) 
which enabled us to deduce, to the accuracy of one //as, the astrometric parameters of a star 
in our Galaxy from observations taken by modern space-born astrometric satellites like Gaia 
(Turon et al. 2005), which are fully consistent with the precepts of General Relativity. 

In this paper we present an astrometric model which contains an extension to the dy- 
namical case, i.e. with the inclusion of the 1/c^ terms, of our previous model which was only 
accurate to 1/c^ (de Felice et al. 2004). We term it as RAM0D3 since it was intended as the 
successor of two previous models (de Felice et al. 1998, 2001). Following the same scheme, 
we shall refer to the model described here as RAM0D4. 

The inclusion of terms of the order of 1/c^ corresponds to an accuracy of 0.1 /tas, at least 
one order of magnitude better than the expected precision of the Gaia measurements. Here 
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we require that the Solar System is isolated and source of a weak gravitational field. These 
two conditions imply that the velocities of the gravitational sources within the Solar System 
are very small compared to the velocity of light, typically of the order of 10 km s~^. Under 
these conditions we select a coordinate system [x\x'^ = ci)(j=i^2,3) such that the background 
geometry has a post-Minkowskian form: 

9a(3 = Vap + ha(3 + (/i^) (1) 

where riad is the Minkowski metric and hai3 are small perturbations which describe effects 
generated by the bodies of the Solar System. These perturbations arc small in the sense 
that \hai3\ -C 1, their spatial variations are of the order of \haf3\ while their time variations 
are of the order of \hap\/c; here and in what follows Greek indeces run from to 3. Clearly 
metric form (1) is preserved under coordinate gauge transformations of the order of h. The 
ha/3^s contain terms of the order of at least 1/c^, hence we shall keep our approximation to 
first order in h. To the order of the time dependence of the background metric cannot 
be ignored therefore the time-like vector field rj — do tangent to the coordinate time axis, 
will not in general be a Killing field (namely an isometry for the space-time) unless one 
moves to far distances from the Solar System where the metric tends to be Minkowski's. 
The components of the vector field rj are f]^ = S^, Vti = 9o^l hence we easily deduce that the 
congruence Crj, namely the family of curves having the vector field r] as tangent field^, will 
have a non zero vorticity. Prom its definition 

U;a/3^P{v)a''P{v)p'^[pVa], (2) 

where P{r])'^ — S'^ + rjxT)^ is the operator which projects orthogonally to rj, Vp the covariant 
p— derivative relative to the given metric and square brackets mean antisymmetrization, we 
have that the only non vanishing components are, to the lowest order: 

Uij = d\jhi]Q O (l/c^) , i, j = 1, 2, 3. (3) 

We remind that /loi and djhQi are ~ 0(l/c^). Condition (3) implies that the surfaces 
t — constant are not orthogonal to the integral curves of r/ at least in and nearby the 
Solar System (de Felice et al. 2004). Nevertheless, because the time lines with tangent field 
r) are asymptotically Killing and vorticity free then the slices S{t) : t — constant allow for 
a non ambiguous 3 + 1 splitting of spacetime with a coordinate representation such that 
space-like coordinates are fixed within each slice. We fix the origin of the space-like coor- 
dinates at the barycenter of the Solar System and assume that the spatial coordinate axes 



These curves are also termed integral curves of the vector field rj. 
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point to distant sources chosen in such a way to assure that the system is kinematically non 

rotating. Adopting the lAU prescriptions (lAU 2000), this system is termed Barycentric 
Celestial Reference System (BCRS) and will be our main system since all coordinate ten- 
sorial components will be relative to it. The non stationarity of the spacetime makes the 
construction of a relativistic astrometric model much less straightforward than for the static 
case considered in de Felice et al. (2004). 

In section 2 these complications are handled to define the geometrical environment for 
light propagation and data analysis. The spacetime metric is not diagonal since terms as 
hoi are different than zero; they arc mainly generated by the velocity of the metric sources 
relative to the given BCRS. Then the gravitational potential at each point of the light 
trajectory depends on the sources in the Solar System at the appropriate retarded position 
as specified in section 3. This will have consequences when fixing the observables and the 
boundary conditions for the differential equations of the light rays in section 4. Section 5 
shows how this model was numerically tested and presents the results of these tests. 

In what follows we shall use geometrized units such that c — 1 — G, G being the 
gravitational constant; with these units a mass Al, expressed in kilograms, has the dimension 
of a length according to the relation M — GM/c^; similarly the time coordinate x° = ct 
will be simply written as t and the spatial velocities are in units of c. For sake of clarity, 
the velocity of light c will appear explicitly only when we specify the order of magnitude of 
terms under discussion. 



2. The light trajectory 

We require that spacetime admits a family of hypersurfaces S{t) with t — constant so 
that the spatial coordinates {x*} are fixed on each of them. As said, these surfaces are 
constrained by the condition of being asymptotically orthogonal to the time direction r) 
which will also be asymptotically Killing. In the nearby of, and of course within the Solar 
System, the spatial coordinates arc not constant along the normals to these hypersurfaces; 
along them, in fact, they vary according to the shift law 5a;* ~ hoidt. To the order of 1/c'^, 
all terms proportional to /loi are in general not zero and cannot be made vanish in a gauge 
invariant way. Let us term u a vector field parallel to 77 and tangent to the family of 
timelike curves 7, say, along which the spatial coordinates are constant; furthermore let & be 
the parameter on these curves which makes the vector field u unitary, namely u°'Ua = — 1. 
Clearly along each integral curve of u carrying the coordinates x', the parameter d[x'^,t) is 



- 5 - 



function of the coordinate time t only. By definition, the vector field u has components: 



e'^5" 



" " da " " "° (4) 

where = (dt/da) = {—goo)~^^'^- The integral curve of u through each spacetime point, 
identifies a local barycentric observer since, as stated, this observer is at rest in the BCRS. 
Consider now a null geodesic T stemming from a distant star at the emission event : 
{xl, t*) and ending at the observation event Pq : (xq, to)- Our aim is to find, from appropriate 
observational data collected at Pq, the coordinate position of the star and its proper motion. 
Let k be the vector field tangent to T and satisfying the following relations: 

kak'^ = (5) 

— + r%,k^k'' = (6) 

where (6) is the geodesic equation, A is a real parameter on T and F"^^^ are the connection 
coefficients of metric (1) given by: 

= ^V"' {di^hp, + d,hp^ - dM + O (h^) . (7) 

Condition (3) holds for the family of the integral curves of u as well, hence Cu does not 
admit a global family of orthogonal hypersurfaces although such a surface always exists 
locally, namely in a small neighborhood of each point. This means that it is not possible 
to define a rest-space of the barycentric observer which covers the entire spacetime. The 
slice S{t()) for example is not the rest-space of the observer u at the observation point Pq- 
In fact, neighboring events belonging to the local rest-space of u at Pq and having spatial 
coordinates x^ = Xq + 5a;*, will have a coordinate time given by 

t ^to + e^'^hojSx^. (8) 

The rest-space of the local barycentric observer u is locally identified by the operator 

P{u)°'i3 = 5^ + u°'Uf3 which projects orthogonally to u. We define the instantaneous line 
of sight of u as the local spatial direction of light propagation; at each point of the photon 
path this direction is found projecting the tangent vector k in the rest-space of u namely 

r = P{uYpkf^ = A;'* + u'^{upk^). (9) 

In this way we define a vector field £ all along the null curve T (see figure 1); the knowledge 
of coupled with that of u"', allows one to determine k"' and therefore to reconstruct the 
null trajectory followed by the photon from the star to the sateUite. It is convenient to 
parameterize the null curve T with the quantity a{\) = &{x{X),t) which marks the proper 
time of the barycentric observer u which the light trajectory crosses at each t. 
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Fig. 1. — The light trajectory, identified by the four-vector fc, propagates in the space-time 
until it is intercepted by the Gaia-like satellite at time Iq. At each point on its trajectory 
the light signal strikes the locally barycentric observer u who identifies in its instantaneous 
rest-space (dotted area) the local line-of-sight £. The surfaces S(t) : t = const do not in 
general coincide with the local rest-space of u. 



If we set: 



UPkn 



we find from (9) that = 1- If we define similarly: 

k"" dx° 



A;" 



u^kp da 



(10) 



(11) 



then the equations which will solve the problem of finding the coordinate position of the star 
at the emission, are fully determined if we set along T 



dx" 
d^ 



- u'- 



(12) 



Prom (6), (7) and (10) the differential equations of the local line of sight read: 

d/" 1 - 1 - - 



+ V 



{d/shxo + dohxi3 - dxhop) 



+ t"[dohxo-^dxhoo 



0. 



(13) 



Let us recall an important property of the vector field £. ^Prom (9) and (4), we deduce 
that Iq — P{u)opkf^ — 0; this imphes 



(^ = ho,t + 0{l/c'), 
hence, to the order of 1/c^ equation (13) becomes: 



(14) 



Afa 1 

— + 6;;td,hoo + -S^i'tdihoo 
da Z 



+ S^) -tfdoh. 



13 



+ r)' 



ef i^dihxj - -dxKjj + 
{dihxo + dohxi - dxh 



Oi) 



r^^dxhoo = 



(15) 
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where i, j, k . . . — 1, 2, 3. Property (14) and the requirement that our analysis is to first 
order in h allows us to decouple the time component £° from the spatial components in 
the set of differential equations (15). We finally have to the 1/c^: 

di^ -- 1 

— -efdihoj--dohoo^O (16) 

Afk 1 / 1 

+ ^H%hoo + 1 (dihko + dohki - dkhoi) 

- Idkhoo = 0. (17) 

The boundary conditions needed to solve equations (12), (16) and (17) are the coordinate 
positions of the sateUite with respect to BCRS at the observation time to and the local 
hne of sight direction l"^^ . While xl are supposed to be known, the components Z(o) need to 
be expressed in terms of specific observables which depend from the experimental set-up. 



3. The reteirded time corrections 

Integration of equations (16) and (17) obviously requires that one calculates the metric 
coefficients ha/s all along the integration path. Due to the linear regime, each perturbation 
term can be written as: 

v = E^S (18) 

a 

where the sum is extended to all bodies of the Solar System. To the selected order of 1/c^ and 
confining our attention to mass monopole terms of the gravitating sources, it is appropriate 
to consider a standard solution of Einstein's equations in terms of retarded tensor potential 
(Weinberg 1972; Misner et al. 1973; de Fehce and Clarke 1990), which can be speciahzed 
also as the Lienard-Wiechert ones (Kopeikin and Schafer 1999). In the last case and in 
conventional units the metric coefficients read: 
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where Wj is a spatial potential which describes the dynamical contribution to the background 

geometry by the relative motion of the gravitational sources and to be defined shortly, 
is the distance from the points on the photon trajectory to the baryccnter of the a-th gravity 
source at the appropriate retarded time. The retarded position of the source is fixed by the 
intercept of its world-line with the past light-cone at any point of the photon trajectory. 
We shall now calculate this distance warning the reader that we shall use geometrized units 
{c— I = G) again. 

Let x'^{a(t)) be the spatial coordinates of a general point P on the light trajectory T 
at the coordinate time t and let x*(cr(t)) be the spatial coordinates of a general point on 
the spacetime trajectory of the a-th source, a being the parameter along that curve. The 
metric coefficients at the point P of the light trajectory are determined by the a-th source 
of gravity when the latter was located at a point Q of its trajectory at a value a{t') of its 
parameter where t' is the retarded time: t' — t — \Ax{t')\ where \Ax{t')\ is the modulus of 
the vector of components Ax'(i') = x''{a{t)) — x^{a(t')). We have dropped the suffix (a). Let 
the world-line of the given spacetime source be described by the tangent vector: 

u"{a) = ^ = -{uaii'^) + 5") , (20) 
da 

where v"" is the a-component of the physical spatial velocity of the source in the rest frame 
of the local barycentric observer u whose components it" are given by (4); hence we have: 

^ = -(upu^)v^ (i = l,2,3). (21) 

The distance at the retarded time t' along the world-line of the source is given by = 
—UffAx^ (sec e.g. Kopcikin and Schafcr 1999). Since we require that our model is accurate 
to the order of l/c^ then, recalling that the metric potentials are at least 1/c^, it suffices 
that the retarded distance is determined only up to (1/c). Prom (20) and recaUing that 
Ax^{t') is taken along the generators of a hght-cone, we have 

r^^\Ax{t')\-5ijv\t')Ax^{t'). (22) 

Before proceeding, let us briefly recall the concept of spatial distance. Let r, say, be such a 
distance from P to the a-th source at Q as it would be measured by the barycentric observer 
u in P. Prom the very definition of a spatial distance, the quantity r is the separation 
between the events P and Q the latter being the event of the history of the observer u going 
trough Q which is simultaneous to P with respect to u in P (see figure 2). 
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Fig. 2. — The metric coefficients in P at time t are determined by the gravity source located 
in Q at the retarded time t' . At time t the gravity source will be in Q but the distance r 
which enters the metric coefficients is that between the events P and Q; the latter has the 
same spatial coordinates of Q and Q but is simultaneous to P with respect to the locally 
barycentric observer u in P. Here we confuse Q with Q the error being of the order of 
In the figure the suffix (a) of the bodies are dropped. Q is the position on S{t) of the planet 
which is moving along its world line. 



Prom the condition of simultaneity, the coordinates of Q are given by: 



x\a{t')) 

rx^ait)) 

t+ / e2%j(x')dx^'. 
JxHHt')) 



Here the integral is at least of the order of 1/c^ and since it would enter terms of an order 
higher than that we can neglect it and confuse Q with the event Q which is the intersection 
of the world- line of u through Q with the surface S{t). To first order in h and terming r 
the parameter along the shortest space-like path (geodesic) connecting P to Q on S{t), the 
required spatial distance reads: 



/■r(P) 
Jr{Q) 



> I 

Jr{Q) ^ 



(23) 



where = dx^/dr is the unitary space-like tangent to the line of integration. Notice that 
the integral in (23) contains terms of the order -\- 0(l/c^)] hence, on the basis of what 
we said before, we can retain only the Euclidean part, namely: 



/ v^;;eFdr+o(i/c2) 

Jt{Qi) 



= |Aa;(OI + C>(l/c2). 
Here we have made the identifications: 

x\Tm 

x\t{P)) 



(24) 



x\a{i!)) 
x\a{t)). 



From (24) and (22) we can write the metric coefficients (19) as 



(a) _ 




00 ~ 




(a) _ 


^2GA<('^) 








o...(«) 



(1 + v(") • n(«)) + O 



(25) 
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where n*^") = ^i") jr^")^ all quantities here being calculated at the retarded time t' . Having 
defined as retarded distance the quantity r-n, we shall term reduced distance the quantity r 
as in (24). 

It is convenient, for computational purposes, to express the distance r in terms of 
quantities defined all at the same time t of metric determination. To this purpose few more 
steps are needed for the correct identification of the retarded time corrections. Integrating 
(21) from d{t') to a{t) which correspond respectively to the points Q and Q as shown in 
figure 2 



x\a{t')) ^ x\a{t)) + / {upu'^)v\a)da--- 

Ja{t') 



(26) 



the integral being taken along the world-line of the gravity source. From this it follows: 



X 



{a{t)) - x'{a{t')) = x'{a{t)) - x'{a{t)) - / {upu^)v'{a) da + ■ ■ ■ . 

Ja{t') 



(27) 



Notice that a;*(a(t)) are the spatial coordinates of the gravity source at the time t in the 
BCRS; these are supposed to be known; also the components of its spatial velocity v'^{a{t)) 
with respect to the local barycentric observer at time t are assumed to be known. Prom (23) 
and (27) we then have: 



r= {6,, 



x'ia{t)) - x'{a{t)) - / {upu^)v'{a) da ■ ■ 

pa(t) 

x^{a{t)) - x^{a{t)) - / {u0u'^)V{a) da + ■ ■ 

Ja(t') 



1/2 



+ 



Setting 

and recalling that 



Ax'{t) ^x\a{t))-x'{a{t)). 



(28) 

(29) 

(30) 



-(m"w«) = (1 - 52)-i/2 

V being the (modulus of the) spatial velocity of the gravity source with respect to the local 
barycentric observer, equation (28) writes as: 







Ax'{t) + 


/ V 




'a{t') V 






Ax^{t) + 






Jait') \ 



r.2 1 



X 



1/2 



(31) 
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which can be written as follows: 



a{t) 



■?/ da 



a(t') 



da 



1/2 



where at any t 



(32) 



(33) 



For the bodies of the Solar System the quantities j^j^^}-^v^ {a) da are always much smaller 
then ro; even in the extreme case of a light ray skimming the surface of a planet the integral 
goes to zero while rg remains of the order of the body's radius. Thus the following Taylor 
expansion is justified. Remembering that (1 + x)^/^ ~ 1 + (l/2)a; — (l/8)x^ + • • then to 
the given order we have from equation (32): 











J a 



v^{a)da + -Sij 



ait) 



ait'] 



8r^ 



a{t') 



v^{a) da + -Sij 



a(t) 



vHa) da 



a{t) 



vHa) da 



ait'] 



a{t') 



v\a) da I v^{a) da 

a{t') 



+ O 



(34) 



Evidently, the 1/c^ contribution of the relative velocity of the gravitating sources to the 
retarded time correction can be neglected. The integrations in (34) contain the unknown 
r implicitly in a{t') and this makes the calculations rather cumbersome, however order of 
magnitude considerations allow one to identify the proper-time a{t) with the coordinate time 
t so achieving a considerable simplification. From (20), in fact, we have 



dt 



da ^ ' 



(35) 



where 7 = —UcrU'^ = (1 — i;^) ^^'^ is the Lorentz factor of the gravity source relative to the 
local barycentric observer. Recalling (1) and (4) we easily get: 



u 



1 + 2^00 + o 



o 



the latter relation arising from the condition VaU"" — 0. Thus equation (35) becomes: 

— = 1 + -/ioo + ^v^ + 
da 11 



(36) 



(37) 
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which leads to 



da^dt + 



(38) 



A similar argument applies to the integration variable a entering the master equations (16) 
and (17) and so, to the required order, we can calculate the individual terms in (16) and 
(17) and in (34) confusing the barycentric and the gravity source proper-times with the 
coordinate time t. Prom the above considerations and setting Ax*(t, f) = x^{t) — x^{t') ~ 
Jl v\t) dt + 0{l/(?) with x\t) = x'{a{t)), equation (34) becomes: 



r = ro < 1 + ^ 



5ijAx'{t)Aff{t, t') + -SijAx'it, t')Aff{t, t') 



8rl 



6ijAx'{t)Ax^{t,t') + l-6ijAx'{t,t')A£'{t,t') 



o\-\. 



(39) 



It is clear that, being v^(f) a bounded function, the integrals in (39) remain finite 
when ro — oo, hence in that limit we have that r — Tq as expected. Let us now define 
ro = {Aa;*(t)} and f = {A5*(t, t')} with i = 1, 2, 3; since the orbits of the spacetime sources 
are bounded then f is finite while ro can grow to infinity, then we retain only terms of the 
order of (f/rg)^. Equation (39) can then be cast in vectorial form: 



r = ro 



1 + — ( ro ■ r + -r • r 



1, 
-1 
2 



^0 



«4 



(40) 



It should be stressed here that (40) is an implicit equation since r is function of r in its turn. 



4. The boundary conditions 

The differential equations (16) and (17) are of the general form: 

d£° 

— = T\dph{x,y,z,t),l\a{x))) (41) 

where are real, non singular, smooth functions of their arguments. A general solution of 
(41) is 

r=r(<7,£7o)) (42) 

where Z^q^ are the components of the vector i at the observation and represent the boundary 
values that we need to fix in order to integrate (41). These boundary conditions can only 
be expressed in terms of the satellite observations. In the case of Gaia the observables are 
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the angles that the incoming hght ray forms with the spatial axes of the satellite attitude 
triad. The latter is a set of three orthonormal space-like vectors which are comoving with the 
satellite and define its rest frame; their coordinate components are denoted as {£^?'}a=i,2,3 
and satisfy the condition E^E^^ = 5^^. The observables are given by: 

^"^^("-^)""^= (PK)./.W)V^ ' (43) 
P{u')af3 is the operator which projects into the satellite's rest-frame, namely: 

P{u')ap = gap + (44) 
and u' is the time-like and unitary satellite four velocity. The latter is given by 

u' = Ts{dt + (5id, + (52dy + (5^d,), (45) 

where ^^'s are the coordinate basis vectors relative to the barycentric celestial reference 
system, (3i are the coordinate components of the satellite three-velocity with respect to the 
barycenter of the Solar System. Here the subscripts refer to contravariant components not 
to confuse them with power indeces. It is important to stress that the modulus of /3j is not 
the physical velocity of the satellite with respect to the local barycentric observer w, nor the 
j3i are its physical components. The physical velocity instead, namely the quantity which 
would be measured, is given by the modulus of the four vector with components: 

i/" = 1 {u"^ - 7m") (46) 
7 

where 7 is the instantaneous sateUite's Lorentz factor. Form (46) it follows: 

= = 0(l/c^) 

= (1 + [7)^ + 0(1/0^) (47) 

where 2U = /iqo Z?^ = /3f hence i/* coincides with /3j only to the order of Finally 
we recall that u'°'u'^ — —1 and this relation fixes the time factor of u' as = 1 + |(/ioo + /?^)■ 
From (12) and (45), equation (43) can be written more explicitely as: 

Piu%piif,^ + u-)E/ (Ifo) - ^^)Eap 

all terms being calculated at the observation time. We easily see that all quantities contained 
in (48) are known except which obviously are the unknown boundary conditions, as 
stated. In order to deduce them correctly one needs the components of the satellite attitude 



-16- 



triad. They have been deduced in exphcit form in Bini et al. (2003) (see equations (4.17) 
and Appendix B of that paper) . Here we shall cast them in a form easier to handle in direct 
calculations. 

The satellite is expected to rotate with an angular velocity LOr about an axis (its x-axis) 
which forms an angle a with respect to the instantaneous local direction to the Sun. The 
spin axis then precesses^ about the direction to the Sun with an angular velocity lOp. The 
attitude triad will then depend on this parameter specification as follows (Bini et al. 2003): 

= C0SQ;Aj(bs) — sinci;cos(tc;pt)A|(bs) — sinQ;cos(a;pt)Ag(bs) (49) 

= — sinQ;sin(cJrt)A"(bs) + 

+[cos(a;rt) cos(c(jpt) — sin(c<Jrt) sin(cL;pt) cosQ;]A?(bs) (50) 
+ [cos(c<Jrt) sin(c<jpt) + sin(a;rt) cos(a;pt) cosa]Ag(bs) 

= — sin a cos(u;r^)Aj (bs) 

— [sm{urt) cos{ujpt) + cos{urt) sin(ci;pt) cosQ;]A|(bs) (51) 

+ [— Sm{u!rt) Sm{u!pt) + COs{u!rt) COs{u!pt) COSQ;]A|(bs) 

where a takes values as {t, x, y, z) and {Aa(bs)} is the Lorentz boosted triad adapted to the 
satellite whose components are given in explicit form in the appendix A. 

4.1. The spatial velocity of the Gaia satellite 

To make the boundary conditions complete we need an operational definition of the 
satellite's spatial physical velocity The sateUite trajectory will be close to that of the 
outer Lagrangian point L2 relative to the Earth-Sun system. The world lines of L2 and Gaia 
never intersect because Gaia moves around L2 in a halo type orbit. Hence to define Gala's 
spatial velocity with respect to L2 we have first to fix a coordinate frame comoving with L2. 
Denoting the coordinates with respect to L2 with a bar, we define: 

i — t — G ^ L 

hence the four-velocity of Gaia in the coordinate frame of L2 reads: 



^In the case of Gaia the sateUite wiU make one turn every ~ 6 hours with a precession period of ~ 70 
days and a precession angle a of ~ 45°. 
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where 

dxi 
^° dt 

are the coordinate components of Gaia's spatial velocity with respect to L2. 

Let us identify what are the unknowns and what is known. Our task is to express the 
formers in terms of the latters. The unknowns are the components of the spatial velocity 
of Gaia with respect to the local barycentric observer, namely z/*. The known quantities 
are: (i) the components of the coordinate spatial velocity of L2 with respect to the local 
barycentric observer namely P}^^ = dx\^/dt] (ii) the components of the spatial velocity of 
Gaia with respect to L2, namely Pg', (Hi) the metric coefficients at the position of Gaia. 

Let us express t/* in terms of the known quantities. From (46) it follows that 

= (^1 + ^hoo + 0(l/c^)^ (52) 

hence, recalling that /3* = dx^/dt, we have 

dxV, dx\ _. . ^ ^ 

namely 

= /9g (1 + Ihoo + 0(l/c^)) + (54) 
In (52) and (54) the gravitational potential /iqo is calculated at the position of the satellite. 



5. Testing RAMOD4 

The most efficient way to test RAM0D4 is to compare the results of the integration 
of its equations of motion with those of RAM0D3. We expect that the differences between 
the two models are of the order of 5 ■ 10^^^ rad i.e. ~ lO/ias (see section 7 of de Felice 
et al. 2004). Obviously the practical implementation of these equations requires the precise 
specification of some of its ingredients. First of all an explicit form of the metric (25) should 
be chosen; in this context and to the prescribed order of 1/c^ we chose, as stated early, a 
solution of Einstein's equations based on the Lienard-Wiechert tensor potentials in terms of 
which, we recall, the post-Minkowskian metric perurbations take the form 

= E^(i--^-)+o(i) 
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= Z^-^M^^i= -2^00- + 0(^1 (55) 

where the potential wj"'* due to the relative velocities of the bodies is equal to 

The coordinate form for the master equations obtained by using the explicit metric com- 
ponents (55) in (17) is reported in appendix B. The reduced distance r*^"^ which enters 
equations (55) will be treated separately in subsection 5.1. 

The boundary conditions which are needed to integrate (17) can be specified only after 
the tetrad describing the motion of the observer is given explicitely. Our goal is to compare 
this new model with RAM0D3, hence a natural choice is a tetrad that makes the measure- 
ments compatible with those described in de Felice et al. (2004), that is, a phase- locked 
tetrad associated to an observer moving on a circular orbit around the barycenter of the 
Solar System on the plane z — Q. 

The four- velocity of this observer writes 

u'- = e^' [5- -uj{y,- y^) +u;{x,- x^) 5!^] , (56) 

where u is the angular velocity of the observer, Xq and yQ are the barycentric position of 
the Sun, and e^' is a normalization factor required by the condition u'°'u'^ — —1. Since in 
RAM0D4 the only non- vanishing terms of the metric are g\a\a and ga, this factor is 

e"^' = {- [goo - (ys - Ve) 9ox + {x^ - Xq) goy + tJ' {y^ - y^f g^^ + u'^ {xs - x^f gyy] ] 

The explicit expression of the components of the tetrad spatial axes can be deduced by 
the components of the tetrad {Aa(bs)} given in eqs.(A3-A19) setting = and substituting 
the velocity components of the observer with the following expressions 

/3i = -uj {y^ - yo) e^' 

(32 = u{x,-xe)e^' (57) 
= 

which describe, as stated, a spatially circular orbit. 
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From (57) and following the conventions in Bini et al. (2003), the tetrad component 
A(bs)j become 



AO 



i(bs) = i '^hoo + ^^^R^c'^'^' + 1) [-a; {y^ - y©) e'^' cos (j)s + uj {x^ - Xq) e^' sin 0^] + 
+ (/iox cos 0s + hoy sin 0^) , 



where R is the radius of the circular orbit of the observer. Since (xg — Xq) — it! cos 0s and 
iUs — He) — -R sin 0s, the above expression simplifies as 

(31 \ UJQ^' 

-hoo + -u'^R^Q^'''' + 1 j [- (Xs - Xq) (ys - Vq) + (a^s - Xq) (ys - Z/©)] + 

[/iQa: (Xs - Xq) + /ip^ (^g - Vq)] 

R 

= [hox {x^ - Xq) + hoy (ys - y©)] ■ 
Similarly, the other components of Aj(bs) turn out to be 



Af(bs) 




R \ 


2 ) 


Aj(bs) 


1 

'r 


{Xs — Xq) — 


2^00 


Af(bs) 


= 0, 







(58) 



while the remaining axes are 



A^ 



2(bs) 
2(bs) 
A|(bs) 
A|(bs) 



= -ic^i? (ys - t/o) C^^' - (1 - ^) 

= lo.2^(xs-Xe)e2'^' + i2^(l-^) 



(59) 



0, 



and 



A°(bs) 
A|(bs) 

A3(bs) 

A|(bs) 



— 

= 

= 

= 1 - 



(60) 



ftoo 
2 ■ 



Finally, given these tetrad components, one can fix the boundary conditions t^^-^ needed 
for the integration of the equations of motion; inverting eq.(48) and to the 1/c^ order, they 
are given by: 

'd 



•^(0) ~ ~ ^ — x,y,z 
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where 



D 



-62/3 ( 1 + 4t/ + -/3^ j + (^1 + 2[/ + -/? 
I eJl + tZ + i/^^l -/^fl + S^+i/? 



a; 
R 

-63 (1 + t/) . 



-e2(l + [/+l/32'j+/5('l + 3C/+^/32 



a; 



- |ei (1 + [/) 



5.1. Implementing the reduced distance 



Let us recall that all the metric coefficients contain the reduced distance of the per- 
turbing bodies, so its practical implementation into the formulae for the tests requires to 
make explicit assumptions on the ephemeris. This was not needed in RAM0D3 since the 
metric is stationary in that model. Moreover, in comparing RAM0D3 to RAM0D4, we are 
at liberty to choose the ephemeris of the Solar System bodies. For sake of simplicity we 
consider here that the perturbing bodies move along circular orbits around the barycenter 
of the Sun- Jupiter system; in this case an approximate analytic formula for the retarded 
distance can be given as explained in the following. 

From the numerical point of view, an error Ar on the reduced distance propagates to 
an error A/iqo on the metric determination 



Ah 



00 ^ 



dh 



00 



dr 



Ar 

Ar ~ hoo — • 



The model should be accurate to the {v/cf order^, so all possible sources of numerical error 
must keep this accuracy. We know that /iqo ~ ('^/c)^, therefore, requiring A/iqo < (v/c)^ 
implies (see eq.(l)) 



Ar V 
— < - ~ 10 

r c 



for the diagonal metric coefficients. Therefore we require that the reduced distance r to the 
a-th planet be known with a relative error of 



^^0.1^, 
r c 



(61) 



•^In this section, where expHcit comparisons with non-geometrized quantities is required, we will use the 
{v/c)"" notation for the order of accuracy instead of the 1/c" one used so far. 
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Prom the last consideration it follows that to keep the numerical accuracy of the tests 
to the order of {v/cf inside the Solar System, it is sufficient to retain only the following 
terms in the expression of the retarded distance (40), i.e. 

r = ro + — (ro • f) . (62) 

To prove this statement let us start taking planets moving on circular orbits around the 
barycenter of the Solar System. Since a — t-\-0 {{v/cY), we have 

x{t) = i?cos(ao + cut) 

y(t) = Rsm{ao+ut) (63) 
z{t) = 



and consequently 



x{t — r) = x{t) cos{ujr) + y(t) sin(c<jr) 

y(t — r) = y(t) cos{ur) — x{t) sm{ur) (64) 
z{t-r) = 0. 



Recalling that 

ro ^ {S,jAx'{t)Ax^{t)Y/' = ^{x{t) - x{t))" + (vit) - mf + ^ - mf (65) 

ro = {x{t) - x{t),y{t) - y{t), z{t) - z{t)} (66) 

f = {x(t) - x{t - r),y{t) - y{t - r), z{t) - z{t - r)} (67) 

and substituting eqs.(64) to (67) into (62) and after some algebra, we obtain the following 
expression (in geometrized units) 

r = ro H [{xx + yy) — {xx + yy) cos(ci;r) — {xy — yx) sin(ci;r) — B? \i — cos (cur)]] , (68) 

^0 

where, for the sake of simplicity, we put x{t) = x and x{t) = x. 

The previous expression can be solved numerically in principle to any degree of accuracy, 
but it cannot be solved analytically because equation (68) is transcendent, like, for example, 
Kepler's equation. 

However, for all the planets in the Solar System, a first-order analytical expression for r 
is sufficient to satisfy the requirements set by (61) about the reduced distance, and therefore 
to have the metric coefficients approximated to (v/c)^. With no loss of generahty, we can 
take a planet whose orbital radius is R, a series of events P = {x, y,z,t), and assume that: 
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1. the spatial coordinates of P all belong to the x-axis (that is P = {x, 0, 0,t)); 

2. the orbit of the planet is equatorial, i.e. z{t) = V t; 

3. the planet has coordinates x(t) — R, y{t) — 0. 

With these hypothesis tq = x{i) — R, and we expect that r — Tq when r/c = At — kT 
{k = 1,2, . . . ,n) where T is the orbital period of the planet. 

Once the series of events P is generated then putting x = R + cAt and using At as the 
independent variable in the range (0, 5T] for each planet of the Solar System, we calculate 
for each P: 

1. the "exact" reduced distance r by solving numerically equation (40); 

2. its zero-order approximation, i.e. ro — {R + cAt) — R; 

3. its first-order approximation obtained from (68), namely 

— roH — {(xx + yy) — ixx + yy) cos(a;ro) — ixy — yx) sin(a;ro) — R?[i — cos (a;ro)]| ; 
ro 

(69) 

4. the relative error for tq, Atq/t — {r — ro)/r; 

5. the relative error for ri, Ari/r — {r — r-ij/r. 

The plots in figure 3 show that, as expected, ro is not a good approximation for r in the 
{v/cf unless P is at least ~ "icT far from the planet, while ri is a good approximation since 
it is always Avi/r < {v/cf. 
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Level approximation for Mercury 



Level 1 approximation for Mercury 



Fig. 3. — Plot of Aro/r (left) and Ari/r (right). The solid line is 0.1 • (t'mercury/c) , i.e. the 



floor for dr/r according to (61), in the left plot, and 



^mercury 



/c)^ in the right one. The 



dotted lines indicate the relative errors. The two plots refer to the case of Mercury and are 
representative of all of the other planets, which follow the same trend. It can be seen that 
the zero-order approximation (ro) is not sufficient to get the required numerical accuracy of 
the reduced distance all along the integration path, while the first order one (ri) keeps the 
differences always well under the {v/cf level. 
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5.2. Test design and results 

It can be easily verified that if we put the velocities of the perturbing bodies to zero, the 
formulae for RAM0D4 reduce to those of RAM0D3. Therefore the first thing checked was 
that the new code produced the same results as the orig inal RAM0D3 code when v(") = 0. 
This cannot be properly called a test, nevertheless it is a basic sanity check to vahdate the 
correct implementation of the model in the computer code. 

After this check, we started the real test phase having in mind that the proper design of 
the tests should highhght the effect of the non-stationarity of the space-time, which marks the 
main difference between RAM0D3 and RAM0D4, as thoroughly discussed in the previous 
sections. 

First of all, the stationary metric of RAM0D3 allowed us to design a straightforward 
test of spherical-symmetry, a property peculiar to such a gravitational field. The same is not 
possible in RAM0D4; therefore identical geometrical configurations adapted for RAM0D4 
should produce differences < 10 //as, which is the intrinsic order of accuracy of RAM0D3. 

Moreover, the non-stationariety of the space-time of RAM0D4 can be conceived as 
being due to the contribution of three terms, namely: 

1. the motion of the planets (i.e. the time dependence of the positions of the bodies r{t)); 

2. the inclusion of the retarded distance, that takes into account the finite propagation 
of gravity (negligible at the order of 1/c^); 

3. the presence of terms to the 1/c^ order which depends explicitly on the velocity of 
the perturbing bodies; these include the off-diagonal terms of the metric hoi and the 
retarded corrections proportional to v^'*'' • n*^") in the diagonal coefficients of the metric. 

The above considerations naturally lead us to design two specific set of tests which are 
described below. The second one, in particular, aims at comparing different versions of 
codes from a numerical point of view. 



-25- 




Fig. 4. — Geometry for the self-consistency test. The observer is ahgned along the Sun- 
Jupiter direction and for each angular distance ■0 four stars at symmetric positions with 
respect to the perturbing bodies are considered. In RAM0D4 the spherical symmetry is 
broken by the orbital motion of the bodies. In the figure only the orbit of Jupiter is high- 
lighted, while that of the Sun cannot be drawn because of its small radius. 
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Table 1. Summary of the results for the self-consistency test. The first column shows the 
angular distance ip from the Sun (in degrees), the second column is the mass ratio of the 
planet and the Sun, Sip is the deflection (identical in all the four quadrants) obtained from 
RAM0D3 measured in arcseconds, the A{5ipi) are the differences in /las between the 
deflection of RAM0D4 and RAM0D3 in the i-th quadrant, according to the schema 

depicted in flgure 4. 



ip (dej 


;) 


mj/rriQ 


dtp 


A(5t^i),A(5^3) 


A(5^2) 


A(5^4) 


~ 0.27(= 




10-'^ 


1.7509921 


< 10-1 


-13.2 


13.2 






10-^ 


1.7507492 




-1.3 


1.3 


5) 




10-5 


1.7507249 




-0.1 


0.1 


)) 




10-6 


1.7507224 




< 10-1 


< 10-1 


1 




10-=^ 


0.4667362 




-1.0 


1.0 






10-^ 


0.4666715 




-0.1 


0.1 


)) 




10-5 


0.4666651 




< 10-1 


< 10-1 


2 




10-=^ 


0.2333501 




-0.2 


0.2 






10-4 


0.2333177 




< 10-1 


< 10-1 


5 




10-3 


0.0932902 




< 10-1 


< 10-1 
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Self-consistency test. Let us remind the original configuration of the spherical symmetry 
test done in RAM0D3, i.e. the self-consistency test (de Felice et al. 2004). In the cited 
paper we imposed the condition that the Sun were the only source of gravity and looked 
for the symmetry of the light deflections with respect to its center. Here we considered the 
similar case where, however, the perturbing field is that of the point-like masses of the Sun 
and Jupiter aligned with the observer; the geometrical configuration relative to this test is 
sketched in figure 4. Here we considered an observer at 1 AU from the Sun and 6.2 AU 
from Jupiter, and a set of stars placed at different angular distances from the Sun. For 
each %l) we have taken four stars symmetrical positioned with respect to the instantaneous 
axis joining the observer to the Sun and Jupiter: the anisotropy among the deflections are 
expected to be more evident on the orbital plane of the two bodies. The results are reported 
in table 1, where we can see that both our predictions were satisfled, i.e.: 

1. dipi and dip^ are equal and coincident with those of RAM0D3 to 0.1 /^as, while 5ip2 7^ 
5ip4^ and different from those of RAM0D3; 

2. the differences between and ^■04 in RAM0D3 and RAM0D4 are of the expected 
order (~ 13.2 /xas for Sun-grazing rays and quickly falling below the 0.1 /xas for angular 
distances ip 5 deg); 

3. the deflection is greater in the case where the Sun trajectory is approaching the photon 
path and smaller in the opposite side, where it is getting farther from it. 

As a final verification for this self-consistency test, we repeated several times the same run 
but with a less massive planet. Obviously, we expected that, as the mass of the planet 
decreases, the results tend to those of the perfectly symmetric case: i.e. as mjup — > 0, the 
velocities of Jupiter and the Sun go to zero as well. Table 1 shows that the defiections are 
perfectly symmetric down to the 0.1 /xas level when mpianet = 10~^mjup. 
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Codes 


Contributions 


r{t) Tret V 


R3a 


Yes No No 


R3b 


Yes Yes No 


R4 


Yes Yes Yes 



Table 2: Schematic representation of the different codes under comparisons. While R4 
stands for the one that implements RAM0D4, R3a and R3b indicate different "flavors" 
of a RAM0D3-like code. Among the specific contributions, r(t) means that we consider the 
equations of motion of RAM0D3, but taking into account the planetary motion in the sense 
described in the main text of the article; r^^t that, when calculating the distances (and the 
velocities) of the perturbing bodies from the photon which enter the metric coefficients, we 
consider the finite gravity speed (i.e. the reduced distance); finally, v means that the equa- 
tions of motion are those of RAM0D4 and that the metric is complete with its terms 
depending explicitly on the velocity of the perturbing bodies. In the original RAM0D3 code, 
which strictly implements the model and its hypotheses, none of these contributions 
are taken into account, while all of them are present in R4 at its full extent. 
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Disentangling the non-stationary space-time effects. To this purpose we have checked 
RAM0D4 against different (somewhat hybrid) versions of the RAMOD framework, where, 
according to table 2, they wiU be caUed in the following: 

1. R4: the code implementing the full RAM0D4 model; 

2. R3b: the R4 code without the 1/c^ terms depending explicitly on the velocity of 
the perturbing bodies. In this case the master equations are reconducted to that of 
RAM0D3, whereas the reduced distances are still fully considered and the planets are 
moving along their orbits; 

3. R3a: the R4 code without both the velocity-dependent terms and the reduced dis- 
tances. In this case we consider, for each step of integration, a different position of 
each perturbing body. These positions are those for circular orbits at time ti of the 
i-th step of integration. 

The importance of these tests, as said above, is numerical rather than physical since, e.g., 
the exclusion of the reduced distance cannot be acceptable in a strictly physical sense. It can 
be important, however, if one is interested in the practical implementation of a particular 
code: in this case what one should care about is the efficiency and the speed of the code 
given the required level of accuracy. A similar comparison was done by Klioner and Peip 
(2003) for the model in Khoner (2003), and in this sense the R3a and R3b cases resembles 
cases Pi and of that paper, respectively. 

In these tests we considered only the presence of the Sun and Jupiter, since all the 
planets can be added with no loss of generality in the same way as Jupiter. Two cases have 
been taken into account, as sketched in figure 5: (i) when the observer lies between the Sun 
and Jupiter, or (ii) when the Sun is between the observer and Jupiter. The first case allows 
to evaluate the specific contribution of the Sun and of Jupiter separately, putting them in 
opposite direction with respect to the observer, while the second is more sensitive in order to 
investigate at which level of accuracy the model feels the non-linearity of the superimposition 
of the two-body gravitational field. 
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Case (i) 



V|/ (0,180°) 



Sun 



Obs. 



\|/ (-180°,0) 



Jupiter 



Case (ii) 



Jupiter 



Sun 

H® — 



Obs. 



\\f (0,180°) 



\|/(-180°,0) 



Fig. 5. — Sketch of the geometrical configurations for the tests using the codes R4, R3b and 
R3a. 
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For both of the cases the results indicate that: 

1. the differences between R3a and R4 are ~ 4.5 /xas at most, meaning that the influence 
of the fixed planetary positions are of < 9 yuas. In particular, the greatest effect is for 
photons which graze Jupiter (~ 4.5 /xas) instead of the Sun (< 0.3 /ias); 

2. the differences between R3b and R4 are tipically less than those between R3a and 
R4, but their ratios vary according to the geometrical configuration of the perturbing 
bodies and the observer, going from, e.g., one order of magnitude to about 1. Once 
again the greatest difference occurs in the case of photons grazing Jupiter, as it should 
be expected since its velocity is far bigger than that of the Sun, and it amounts to 
~ 1.8 /xas while the Sun contributes with ~ 10~^ /xas at most. 

This means that one should be careful at choosing the code integrating the photons' geodesies. 
Optimizations are possible, but they must be tuned with the geometrical configuration of 
the perturbing bodies. 
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Table 3. Excerpts from the results of the code-comparing test. This table is for the (i) 
case, in which the observer lies between the Sun and Jupiter. The first column gives the 
angle that the photon's incoming direction forms with the x axis counted counter-clockwise 
in the range [—180°, 180°) as sketched in figure 5. The next three columns contain the 
deflections given by the three codes (in arcseconds) and the fifth and sixth columns give 
the difference in /las of R3b and R3a with R4 respectively. 







<)!.■ (arcscc) 








(/'as) 


(<)l,-R4 - ''><.'R:ii, ) (/'as) 




R4 


R3b 




R3a 








179.7334367 


1.75073792543 


1.75073781174 


1, 


.75073764027 





.11369 


0.28516 


179.4668735 


0.87535666435 


0.87535660782 


0, 


.87535652209 





.05653 


0.14226 


178.4006205 


0.29176702035 


0.29176700193 


0, 


.29176697336 





.01842 


0.04699 


160.1792868 


0.02330937402 


0.02330937317 


0, 


.02,330937094 





.00086 


0.00308 


120.1507651 


0.00707469758 


0.00707469781 


0, 


.00707469723 


-0 


.00024 


0.00035 


80.1222434 


0.00342350034 


0.00342350079 


0, 


.00342350051 


-0 


.00045 


-0.00018 


40.0937217 


0.00148347654 


0.00148347699 


0, 


.00148347671 


-0 


.00044 


-0.00016 


5.7835603 


0.00018740032 


0.00018739969 


0, 


.00018739782 


0, 


.00062 


0.00250 


0.03912 


0.0025469505 


0.0025471242 


0, 


.00254738571 


-0 


.17370 


-0.43521 


0.01956 


0.00480767427 


0.00480800257 


0, 


.00480849601 


-0 


.32830 


-0.82174 


0.00652 


0.01176284745 


0.01176365151 


0, 


.01176485859 


-0, 


.80406 


-2.01115 


-0.00652 


0.0263647168 


0.02636291438 


0, 


.02636020825 


1, 


.80242 


4.50856 


-0.01956 


0.0062147205 


0.00621429593 


0, 


.00621365819 


0, 


.42457 


1.06231 


-0.03912 


0.00289458627 


0.00289438879 


0, 


.00289409164 


0, 


.19749 


0.49463 


-5.7835603 


0.00018738265 


0.00018738327 


0, 


.00018738515 


-0 


.00062 


-0.00250 


-40.0937217 


0.00148347648 


0.00148347603 


0, 


.00148347631 





.00044 


0.00016 


-80.1222434 


0.00342350043 


0.00342349998 


0, 


.00342350025 


0, 


.00045 


0.00018 


-120.1507651 


0.0070746969 


0.00707469666 


0, 


.00707469725 


0, 


.00024 


-0.00035 


-160.1792868 


0.02330936614 


0.02330936699 


0, 


.02330936922 


-0 


.00085 


-0.00308 


-178.4006205 


0.29176611396 


0.29176613238 


0, 


.29176616095 


-0, 


.01842 


-0.04699 


-179.4668735 


0.87534872947 


0.87534878599 


0, 


.87534887173 


-0 


.05653 


-0.14226 


-179.7334367 


1.75070641203 


1.75070652571 


1, 


.75070669718 


-0 


.11368 


-0.28515 
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Table 4. Excerpts from the results of the code-comparing test. This table is for the (ii) 
case, in which the Sun is between the observer and Jupiter. The table headings are the 

same as in table 3. 



?/) (deg) 






Stp (axcsec) 






{'5V'R4 - iSViRSb) (Mas) 


(<5'/'R4 - '5V'R3a) (Mas) 






R4 


R3b 




R3a 






179.7334367 


1, 


.75097892390 


1.75097901914 


1, 


.75097916265 


-0.09525 


-0.23876 


179.4668735 


0, 


.87548418970 


0.87548423717 


0, 


.87548430900 


-0.04746 


-0.11930 


178.4006205 


0, 


.29181111968 


0.29181113523 


0, 


.29181115918 


-0.01554 


-0.03950 


160.1792868 


0, 


.02331311678 


0.02331311760 


0, 


.02331311946 


-0.00082 


-0.00268 


120.1507651 


0, 


.00707631845 


0.00707631834 


0, 


.00707631882 


0.00011 


-0.00036 


80.1222434 





.00342512778 


0.00342512747 


0, 


.00342512763 


0.00032 


0.00016 


40.0937217 


0, 


.00148624154 


0.00148624114 


0, 


.00148624118 


0.00040 


0.00037 


5.7835603 


0, 


.00020574544 


0.00020574502 


0, 


.00020574502 


0.00043 
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0.00043 
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1, 
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1, 


.75100509907 


0.09560 


0.23963 
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6. Conclusions 

Aim of modern relativistic astromctry is to produce a three dimensional rendering of 
our Galaxy with an accuracy of one fias in the measurements of angles. In order to reach this 
goal one needs an appropriate algorithm to trace back to their emission space-time points 
the light signals which are intercepted by a suitable observational device. 

In the near future, one of such devices will be the astrometric satellite Gaia which will 
orbit the Sun from nearby the Sun-Earth outer Lagrangian point L2. The expected accuracy 
of the satellite observations requires that one should treat the light propagation through the 
Solar System in a general relativistic framework and consider the general relativistic effects 
induced by the bodies of the Solar System up to the order of 

Our strategy was to construct a series of models with increasing generality and com- 
plexity, exploiting each of them as a test-bed for the more advanced one. Our last model 
accurate to 1/c^ is RAM0D4. As previously discussed, we first compared RAM0D4 to its 
less accurate predecessor RAM0D3, which was accurate to With a suitable handhng 

of the items considered in the model, we are able to highlight the relative importance of the 
individual sources of relativistic perturbations and judge about their physical relevance. 

Our first conclusion is that RAM0D4 behaves as expected with respect to RAM0D3, 
that is, the differences between the two models are at the expected level of ~ 10 //as. 
So this work successfully completes the series of models from a theoretical point of view 
and represents a launching platform to deduce with a completely numerical algorithm the 
astrometric parameters of a celestial object from a well-defined set of relativistically measured 
quantities. 

This model, however, is part of a broader project whose completion requires its imple- 
mentation into a feasible and efficient structure for the data reduction of the Gaia mission. 
In this sense, the practical application of this model requires also to investigate the specific 
ways of implementing it with respect to the accuracy goals of the mission. For this reason, 
then, we tried to single out the different contribution which come out together from the 
motion of the perturbing bodies, and evaluate them separately. We found that the contribu- 
tion to the total deflection coming from the velocity-induced terms of the metric is generally 
smaller than that of the reduced distance, but not neghgible at the 0.1 //as level. The former 
one can in fact reach a maximum amount of ~ 1.8 /tas while the latter goes up to ~ 4.5 /tas. 
Both these cases, and the greatest differences with RAM0D3, come out when the photon 
is grazing Jupiter, while the Sun contributes up to less than 0.3 /tas that is one order of 
magnitude lower. 

Given the unprecedented accuracy reachable by Gaia-like missions, several new rela- 



-35- 



tivistic tests related to the light propagation in a non-stationary gravitational field could be 
carried out. From this point of view it is advisable to create a cross-check framework of mod- 
els to validate the relativistic outputs induced by the modeling itself at the microarcsecond 
level. In fact, our Relativistic Astrometric Model naturally confronts itself with the works 
of Kopeikin, Schafer and Mashhoon (Kopeikin and Schafer 1999; Kopeikin and Mashhoon 

2002) and Klioner (2003). While the latter stems from the general formulation of the former 
conforming to its basic tenets, our model responds to an inherently different strategy. As 
far as the astrometric problem is concerned, in the Kopeikin-Schafer and Klioner approach 
the light ray is reconstructed as a sum of terms which allows for a direct evaluation of the 
individual relativistic effects induced by the gravitating bodies that the hght ray is expected 
to encounter on its way. Each relativistic effect enters as a perturbation of the light ray 
trajectory and is treated to a sufficiently high order (not less than to cope with the 
accuracy expected from Gala's observations. Our model instead aims to determine a full 
solution for the light trajectory which naturally includes, in a curved space-time, all the in- 
dividual effects; the latters are somewhat hidden in the covariant formalism of our approach 
and directly contribute to the solutions of our master equations. Evidently any specific effect 
one is interested to explore can be independently deduced from our formalism as a branch 
output, once we adapt the master equation to a required specific case. 

Although our model and Klioner's are of comparable accuracy as it can be deduced 
from individual tests, the overall structure of Klioner's model makes it ready to handle 
the reduction of a large number of observational data in a computational efficient way; 
nevertheless we are confident that our model will soon reach the same versatility and operate 
as an essential tool of verification. However, the general covariant formulation of our model, 
including the complete relativistic treatment of the satellite attitude (Bini and de Felice 
2003; Bini et al. 2003), represents a well defined framework where any desired advancement 
in the light tracing problem using astrometric data is contemplated. 

This work has been partially supported by the Italian Space Agency (ASI) under con- 
tracts ASI I/R/117/01, and by the Itahan Ministry for Research (MIUR) through the COFIN 
2001 program. Vecchiato acknowledges the support of the National Institute for Astrophysics 
(INAF) and Crosta the Astronomical Observatory of Torino (INAF-OATo). 

A. The attitude frame 

Given the tetrad adapted to the local barycentric observer {u^ Aa} (Bini and de Felice 

2003) , first we identify the spatial direction to the Sun as seen in the local BCRS and at a 
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general point on the satellite's space-time trajectory. Then we construct the triad {Ao(s)} 
where one of those vectors identifies the Sun direction. Clearly the set {u, \a{s)} forms still 
an orthonormal tetrad adapted to the local barycentric observer. 

As second step, we boost the vectors of the triad {Aa(s)} to the satellite rest frame; 
remembering that u' is the vector field tangent to Gala's world-line we have (Jantzen et al. 
1992): 



7 



7 + 1 



(Al) 



a=l,2,3 



The tetrad {AQ(bs) = u', Aa(bs)} represents a CoMRS Sun-locked frame, i.e. a Center-of- 
Mass Reference System (Bastian 2004) comoving with the satellite and with one axis fixed 
toward the Sun at any point of its Lissajous orbit around L2. The relation between the 
components z/" of the spatial four-velocity u and the components appearing in (45) is 
easily established from (45) itself and (46) and reads: 



7 



(A2) 



The explicit expressions of the components of these vectors relative to the BCRS are reported 
in Bini et al. (2003). 



The triad components are given by 

1 



A 



A? 



i(bs) 
i(bs) 

Aj(bs) 



Af(bs) 



A* 



2(bs) 
I{bs) 



A|(bs) 



A" 



2(bs) 



A 



3(bs) 



^'(1 H — + 3U) + cos 9 cos + cos 9 sin 0/io2 + sin 9h, 

cos 61 cos 0(1 -U) + ^/Si^f 

cos 61 sin 0(1 -U)\ \(^'2^ 

sine(l-C/) + ^/33* 

$(1 W) - sin0/ioi + cos0/io2 

z 

-sin0(l-C/) + ^/3i$ 

COS0(l-C/) + ^/?2$ 
z 

6 ( 1 -I- - -I- 3[/) — sin ^ (cos 0/ioi + sin 0/io2 ) + cos ^/ios 
z 



03 



(A3) 
(A4) 
(A5) 
(A6) 
(A7) 
(A8) 
(A9) 
(AlO) 
(All) 
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A|(bs) = -cos^sin^(l - [/) + ^/?ie (A12) 
A|(bs) = -sin^sin0(l-[/) + ^/?2e (A13) 
A|(bs) = cos^(l-^7) + ^/33e (A14) 



where 



^ = cos 6* cos + cos 61 sin 6'/32 + sin 6'/33 (A15) 
$ = - sin 0/^1 + cos 0/^2 (A16) 
© = - sin 6'(cos0/3i + sin 0/32) + cos 6'/33. (A17) 

Here 2C/ = /iqo and the angles ^ and are defined as 

= tan-^ e = tan-i , (A18) 

V^'l + V% 

where Xq, and z'q fix the spatial coordinate position of the Sun relative to the satellite at 
each time t\ clearly in terms of the coordinates of the Sun and the satellite they are defined 
as 

x'q = Xq-Xq y'^ = yQ- I/O z'q = Zq - Zq. (A19) 

See Bini et al. (2003) for details. It is straightforward although tedious to show analytically 
that £q as solution of (48) with triad (49), (50), (51) is unitary as expected. 



B. Explicit form for the master equations 

Adopting the metric given in (25), the differential equations for the spatial components 
t of the null geodesic take the following form: 



dt2 



1 /d,T^'9r(") ar(") 



2 V dt dt 

3 dx^ 



2 dt 



dx^ 



(1 + 



/ da; 
Vdt 



dy 
dt 



dx dr^"^ dy c^r'-"-* dz dr^""^ 
^ ^ \ 

dt dx dt dy dt dz 



dx'' dr^"^ 1 dr^"^ 
+ 



^ ^ dx^ dt ^ dt " / 



j.{a) 



dt dt 2 dx^ 
3 da;^' da; 9 (v^") • n^")) ^dyd (v(") • n^")) ^dzd (v(") ■ n^'')) 



2 dt \ dt 



dx 



dt 



dy 



dt 



dz 



+ 



-38- 




4m(") dr , 

—r-s ^ • (V X V 

r(«) dt ^ 



In the case of circular orbits, the last line of each equation disappears since V x v*^") = 
in this case. Moreover with straightforward calculations, similar to those appeared in 
section 5.1, one can easily find that 

v(") • n(") = [{yx^"^ - cos(a;r(")) + (xx^"^ + yy^''^) sin(a;r("^)] 

and 



[cos(cur('*))(5fa;(") - Sfy^"'^) + sin(cur("))(5fa;(") + S^y^"^)] + 



+(9ir(")) 1^ [-(yx^") - sin(u;r(")) + (xx^") + yy^"^) cos(u;r(«))] - 

where all the coordinates of the planets are at time t, while r^") and v^"^ are at time t' (this 
means that one should replace r^"^ with its approximation of eq.(69). 
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